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Small cosmological signatures from multi-brane models 
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(Dated: July 15, 2008) 

We analyse the signatures of brane inflation models with moduli stabilisation. 
These are hybrid inflation models with a non-trivial field-space metric which can 
induce complex trajectories for the fields during inflation. This in turn could lead 
to observable features on the power spectrum of the CMB fluctuations through de- 
partures from near scale invariance or the presence of isocurvature modes. We look 
specifically at multi-brane models in which the volume modulus also evolves. We 
find that the signatures are highly sensitive to the actual trajectories in field space, 
but their amplitudes are too small to be observable even for future high precision 



CMB experiments. 
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universe gives us the best possibility to test string theory models. Brane- 



y universe 



and modular inflation models 



6| give predictions that are largely 



compatible with the present cosmological data. Therefore, in order to differentiate between 



the various models we have to look at their predictions beyond the usual set of experimentally 
determined parameters (the amplitude of the density perturbations, the spectral index and 
the running of the spectral index). 

With the advent of high precision cosmological Cosmic Microwave Background (CMB) 
and Large Scale Structure (LSS) data the possibility of characterising the initial perturbation 
spectrum beyond these simplest spectral parameters is now a reality. Efforts to reconstruct 
the initial perturbation spectrum with model-independent 'inversion' techniques are already 
under way (see e.g. J7J) and have yielded some tantalizing hints of structure in the spectrum. 
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There is a lack of power on the largest scales 8| and some indications of oscillations on 
intermediate scale [9j. Arguably, the statistical significance of these effects are not sufficiently 
strong to motivate a departure from the simplest phenomenological models of single-field, 
slow roll inflation but an accurate analysis of observable features in theoretically driven 
scenarios is certainly warranted. 

A generic feature of models of inflation within warped compactification mechanisms is 
the presence of many extra degrees of freedom in addition to the field driving inflation. 
Depending on the masses and couplings of the extra fields their presence can induce non- 
trivial trajectories in the field space. This can result in observable effects such as radically 
broken scale invariance, presence of isocurvature (or entropy) perturbations and enhanced 
non-Gaussianity of the perturbations close to points of high acceleration in the trajectories 

0. 

A non-vanishing isocurvature contribution to the total initial perturbations is a particu- 
larly interesting possibility since tight limits on any such contribution will be available when 
CMB polarization measurements becomes precise. The polarization data will also help break 
fundamental degeneracies in our ability to accurately constrain any broken scale invariance 
beyond a simple running of the spectral index. 



In this work we specifically look at multi-brane models |ll|]. These models are very 
promising since they naturally avoid a number of significant fine-tuning problems present 
in the simplest brane inflation models. Firstly they naturally produce a larger number of 
e-folds via the assisted inflation mechanism [12|], and secondly, they offer the possibility that 
the inflaton potential is flattened dynamically. 

In these models, all branes are initially initially in a local minimum and as branes tunnel 
through the potential barrier and annihilate with anti-branes the barrier is lowered dynam- 
ically until it completely disappears. The potential becomes monotonic and very close to 
flat once a sufficient number of branes-anti-brane pairs have annihilated and the remaining 
branes will roll down the potential before annihilating. Quantum fluctuations will cause the 
rolling branes to start from slightly different initial positions and therefore follow different 
trajectories. We study the signatures that these models give which may be detected through 
cosmological observations. 



All models under consideration feature moduli stabilisation for both the shape [13] and 



volume 
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moduli and anti-branes needed to lift the anti-deSitter vacuum generated by 



the stabilising mechanism to a deSitter one during inflation. The bulk contains a warped 
throat, the anti-branes being located at the bottom of it, and the mobile branes roll down 
this throat while inflation takes place. 

Specifically, we follow the evolution of a number of perturbation modes for a set of scalar 
fields (in this case the positions of the mobile branes) coupled to gravity both inside and 
outside the Hubble horizon. The non-trivial shape of the potential and Kahler metric leads 
to a residual evolution of the scalar fields after horizon crossing. We evolve the perturbations 



of the scalars coupled to gravity using the Mukhanov- Sasaki [15|, [16| variables and decom- 
pose the perturbations into adiabatic and entropy ones. We finally compute the separate 
spectra for the adiabatic and entropy perturbations and we find that while the adiabatic 
perturbations are insensitive to the trajectories followed by the many inflaton fields, the 
entropy ones are highly dependent on the trajectory. However, their amplitude is much too 
small to lead to observable features in the measured CMB spectrum. 

The study per formed here applies to a larger class of models, namely the "inflection 



point" models 17J, where the mass of the inflaton(s) is large, except for a very small region 
around an inflection point of the potential seen as a function of one inflaton field at a time 
(keeping all other fields constant). Most of the e-folds are coming from this small region of 
the field space where the inflaton trajectory can be very well approximated by a straight 
line therefore leading to a very small amplitude for the density perturbations. 

II. THE MODEL 

As mentioned in the introduction, the model we study here covers a more general class 
of brane-anti-brane inflation models. The important features are: 



• Stabilisation of the complex structure moduli via fluxes 13] and of the volume modulus 
via non-perturbative effects [14J, ll8] ) 

• The presence of anti-branes that lift the vacuum to deSitter during inflation. 

• A number of mobile branes that will roll down the warped throat and annihilate with 
the anti-branes, ending inflation. 

In these scenarios the dynamics of the inflaton is determined by the brane anti-brane in- 
teraction and by the mechanism responsible for the stabilization of the volume modulus. 



The vacuum energy responsible for inflation is provided by the (warped) anti-brane tension. 
The stabilization of the volume modulus is done by D7 branes spanning the 3 large space 



CY.D7 





FIG. 1: Brane setup and field trajectories. The potential features a local minimum for the volume 
modulus and an inflection point for the inflaton, the brane-anti-brane separation. 



dimensions and wrapping four-cycles of the compact space. Gaugino condensation in the 



stacks of wrapped D7 branes generates a non-trivial 
which in the simplest case takes the form given in 14]. 



superpotential for the volume modulus, 



W = w + Ae 



-iaT 



(1) 



The non-trivial embedding of the D7 branes in the compact space makes the coefficient A 
dependent on the position of the mobile D3 brane 19|. An additional dependence of the 
scalar potential derived from Eq. ([I]) on the position of the mobile brane comes from the non- 
trivial Kahler metric [5j The resulting effects of the moduli stabilization and the uplifting 
of the vacuum from anti-deSitter to deSitter by adding anti-D3 branes give a potential that 
is g enerally too steep for inflation, the value of the slow-roll parameter r\ being r\ = 2/3 
201 ] . However, when taking into account the competing effects of the moduli stabilization 
and anti-brane attraction, the resulting potential for the mobile brane features an inflection 
point [ll|, 



17 



211 ] . The vanishing of the r] parameter at this point makes it possible for 



inflation to work, most of the e-folds coming from a small region of field space. 

Apriori, such models require a large amount of fine-tuning, as the slope of the potential, 
and therefore the e slow-roll parameter, at the inflection point must be very small. However, 



multi-brane models offer the possibility to alleviate this fine-tuning via the dynamical flat- 
tening mechanism for the inflaton potential. In what follows we will study the cosmological 
signatures of multi-brane models going beyond the calculations of amplitude of the density 
perturbations, the spectral index and it running to determine what unique signatures one 
can expect to be detectable in cosmological observations. 



A. The Kahler potential 



The Kahler potential proposed by DeWolfe and Giddings 22j has the following form: 

/C = -31og(r) = -31og(T + T-/(^)), (2) 



where we will assume that for a multi-brane setup the function / decomposes on a sum of 
functions, each one depending on one brane position ip 1 : 



K = -Slog (t + T -J2fi(<p\¥)) 



(3) 



The calculations in this paper use the simplest form of the functions ff. 



fi {<p*, v*) = <pY- 



(4) 



One can justify this assumption by taking into account the fact that most e-folds of inflation 
are produced around the inflection point of the potential and therefore the functions fi can 
be expanded to the lowest order in the brane positions. 

In the general case the Kahler metric derived from the Kahler potential of Eq.([3]) has the 
form: 
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(5) 



9fn 9f„ 1 8 2 f n 



df n dip n 

The inverse can also be calculated in the general case (see appendix I A ip . We want to write 
the kinetic energy in terms of the real and imaginary parts of the fields [2l| . The imaginary 
parts will play no role in our subsequent analysis as they will roll very quickly to their 



6 

respective minima and these minima do not shift during the rolling of the real parts [18]. 
We write the volume modulus as T = a + iu and the brane positions as if 1 = (j) 1 + i\ l which 
allows us to read off the field space metric and it inverse for the a and <ft fields: 

(6) 





B. The potential 

In this section we describe the various contributions to the scalar potential. They are the 
F-term responsible for the stabilization of the volume modulus, the anti-brane contribution 
that uplifts the resulting anti-deSitter minimum to a deSitter one, and the brane-anti-brane 
attractive potential which gives a non-trivial evolution to the inflatons. The F-term is 
calculated with the usual formula: 

V F = e K (K^DJVDbW — 3 |M^| 2 ) , (7) 

where K, is given in Eq. ([3]). Let us take for example a scenario with two mobile branes; the 
Kahler potential is: 

K = -31og(r) = -3 log (T + T-h (f^) - f 2 (frf)) . (8) 

We take the superpotential to be the racetrack one, this way the vacuum with no anti-brane 



present is already Minkowski 181 ]: 



W = w + Ae~ al + Be- bl . (9) 

If we take the minimal form for the functions f\ ((f) 1 , (j) 1 ) and f 2 ((j) 2 ,(j) 2 ) and the general 

expression for the F-term, Eq. (1A12I) . the resulting potential has a very simple form: 

e~ acr 
V F = —y [aA 2 e~ au (3 + aa) + aAw cos (au)] + 



3/' 

,-fecr 



[bB 2 e- bu (3 + ba) + bBw cos (few)] + 



3r 2 

[3 (a + b) + 2aba] cos ((a - b) lu) (10) 

oT 

The minimal form of the functions f\ and f 2 brings an important simplification as the 
dependence of the F-term potential on the positions of the mobile branes comes only through 

r = T + T- (b 1 ^ - 2 02. 



The anti-branes that lift the minimum of the potential to a deSitter one during inflation 
have a contribution of the form: 



v; 
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E; 



(11) 



The coefficients Ei are the warped anti-brane tension. The exponent a = 2 for anti-branes 
located in a highly warped region, i.e. at the bottom of the throat 20|. It is not necessary 



231 ]. but we will consider a single-throat 



that all branes are located inside the same throat 
model for simplicity. 

We also have to add the Newtonian attractive pieces for the branes and anti-branes: 

M,N 



V, 



K': 



N 



V - 

4/ =1 r^ (^ - $) 



l- 



(12) 



We want to trace the evolution of more than one brane, and the above potential is singular 
when an brane collides with an anti-brane. Therefore we want to regulate the Newtonian 
potential such that it reproduces the correct inverse-power-law at large distances, but at the 
same time is regular and cancels the brane-anti-brane tension at zero separation. We follow 
here the method of Ref. 



llj and for a given brane-anti-brane pair we choose 

M,N 



v 1 



N 



■) 



E 



k 



i;i 



i=i,j=i r< 



8 v + tyi - <f? y 



(13) 



where the constant s u is chosen such that when <p l = (f) J , the Newtonian potential exactly 
cancels the contribution to the potential of the anti-brane placed at 0q. H there are more 
than one anti-brane located at 0q, we choose s y to cancel the tension of only one of them. 
To summarise, the choice of Kahler potential and the assumption that we have a single- 
throat model, such that all anti-branes are located at the same point, allows us to write the 
scalar potential for a model with n mobile branes and n anti-branes as: 



v{?rf...<ir) 
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F(a) + nT-J2 
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s* + (d) 1 - 



(14) 



i=1 - . vr ^0) 

The function F (a) contributed by the F-term is responsible for the stabilisation of the 
volume modulus at some value o"o, and, as long as we do not add too many anti-branes to 
destabilise this minimum, we can simply approximate F (a) = \vn 1 (J (a — a ) when studying 
the cosmological signatures of the model. 



III. PROPERTIES OF THE FIELD SPACE METRIC 



We want to get a better understanding of the properties of the target (field) space. In 
order to do so we would like to write the Kahler metric for the real parts of the fields in a 
more convenient coordinate system. 

dS 2 



— [da 2 - 2((Jd(i) i d<j + re 



+ 



WAJJAfrlAJ 
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(da-^dcj) 1 ) 2 + rc 



6 



-dr + re 



We can now redefine the coordinate r as r = u 2 . This brings the metric in the form: 



dS 1 
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2^„,2 



Au du 



+ u< 
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u- 



du 2 + ^2 ' 



t=i 



(15) 



(16) 



which is the metric for the n+1 - dimensional Euclidean AdS space. The original coordinates 
cover only a limited patch of the space, the boundary of it being given by the condition: 

n 



2-E 



>0, 



(17) 



»=i 



which translates into the condition that u > 0. This is the Poincare half-space model of 
AdS, 24(. The expressions for the metric connection coefficients are much simpler in the 



new variable. Up to symmetries in permuting the indices: 



TTU 



1 

u 



-pu _ ^y 



"j 



u 
u 



(18) 

All other connection coefficients, most notably P - fc vanish and the non- vanishing ones are 
now independent of the brane positions. As far as the potential of the effective theory is 
concerned, one has to replace a by the new variable u: 

2 l(2 

l , a = - \ u 

In KKLT type of models the potential has the form Eq. (TT4l) . this change of variables brings 
the potential to the form: 



r = 2a-£^ = M 2 , <r=^[u 2 + J2^)- (19) 
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i=l 
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The volume-stabilising function F preserves the full SO (n) rotational symmetry of the 
brane-position inflatons, and the dependence on the brane positions is removed from the 
overall factor l/u A . Combined with the fact that in terms of the new variables the field- 
space metric is independent of the positions of the branes, the Newtonian brane-anti-brane 
attraction remains the only part of the potential which sets a preferred direction in the field 
space. This will lead to a field trajectory during slow roll that will be a straight line, and 
therefore to very small signatures in the CMB power spectrum. 

IV. FIELD EVOLUTION AND INFLATION 

To study the dynamics of the fields during inflation we start with the action for the scalar 
fields coupled to gravity: 

S = J d 3 xdtV=9 (—Gijd^^ -V+ ^\ . (21) 

We assume that the space-time metric is of the FRW type, 

ds 2 = -dt 2 + a 2 (t) (dr 2 + r 2 dn 2 ) , (22) 

and we decompose the fields into a background, homogeneous part, and a space-dependent 
perturbation. The equation for the homogeneous fields is the geodesic equation in field space 
with a driving force given by the potential gradient and a damping force given by Hubble 

expansion: 

BV ■ ■ 

<p k + 3H(j) k + G ka — + rJ.0V J = 0. (23) 

The Hubble constant and its time derivative are given by: 

H 2 = ^(±G ij <j> i P + v\ , H = -(AnG)G^^. (24) 

We write here the more general equation for inhomogeneous fields, as we will need it later 
when studying the perturbations around the homogeneous background: 

' -A [V=99^d u <f> K ] - G KA ^- + rfjtTd^W = (25) 



H 2 = \ [Gjjg^d^dA 3 + V) (26) 
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V. PERTURBATION EQUATIONS 

In this section we analyze the perturbations around the background, slow-roll evolution. 
First, the background fields and metric are only dependent on time, while the perturba- 
tions are dependent on space as well. We will follow Ref. 25j and parametrise the metric 
perturbations as follows: 



9v» + <5fiV = 
and that of the fields 6 1 as: 



(l + 2A(t,x)) a(t)B ti (t,x) 

a (t) B ti (£, x) a 2 (t) [(1 - 2^ (£, x)) c% + 2E >ij {t, x)] 



(27) 



(f) I ^(f) I (t)+5(f) I (t,x.). (28) 

We now expand the evolution equation fT23|) for the fields (p 1 around the background. We 
will look at each term in the equation separately. 

A. Kinetic Term 



We start with the first term in Eg. (1251) and write the scalars and the metric as in 
Eq. (l27ll28p . The result for the plane wave expanded perturbations (see appendix IA 21 for 
details) is: 
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(29) 



where k = \k\ is the expansion wavenumber. 



B. Potential and Christoffel Terms 



The perturbation of the potential term is the simplest one, since it does not involve the 
space-time metric. We simply take the variation with respect to the fields (f) 1 : 
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(30) 
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Finally, we have to calculate the perturbation of the connection coefficients coming from the 
non-trivial field-space metric. They are: 



1 u 



+ 6(f)) <T 3„ {(f) 1 + 6(f) 1 ) d v ((f) J + 6(f) J ) = Tfjg^d^d^ 
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(31) 



background 

We can now collect the perturbation terms and write the equation for the fluctuations of 
the fields: 
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(32) 



Using the equation for the background fields, the perturbation equation can be rewritten as: 
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VI. METRIC PERTURBATIONS 



In this section we analyse the perturbations of the Einstein equations, namely the pertur- 
bations of the Einstein tensor and stress-energy tensor. We start with the parametrisation 
of the metric perturbations given in Eq. (127|) : 

/ -2 A oS i aB >2 aB j3 \ 

aB tl 2a 2 (-2ip + E <n ) 2a 2 E <12 2a 2 E <13 

aB, 2 2a 2 £ 21 2a 2 (-2</> + £ 22 ) 2a 2 £ 23 

\aB t3 2a 2 E M 2a 2 E t32 2a 2 (-2^ + E )33 ) J 
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(34) 
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The functions A, ■?/>, B, and E are in general dependent upon all coordinates, t, x, y, z. We 
follow here Ref. [15] and perturb the Einstein equations written in the form: 



R v u - -R5" = &<nGT v a . 

fi 2 A* P- 



(35) 



First, the perturbations of the Einstein tensor is: 
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(36) 
(37) 



For the stress-energy tensor the general expression is: 
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(38) 



We now perturb the stress-energy tensor according to our parametrisation of the metric and 
field perturbations we obtain: 
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(39) 



We can use here the fact that the background fields are homogeneous to further simplify the 
expression for 5T®. We obtain: 
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(40) 



For the component Tq the background value is trivial, Tq = and the only term that 
contributes is the one in the next-to-last line in Eg. (1391) : 



5T = Guftdi. 



(41) 
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Therefore the equation for the evolution of the metric perturbations become: 

2(ij> + -A) =8irG-Gu<j> I 6<l> J i (42) 
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(43) 



2 d<f) K r r r <90 x r 

VII. THE EVOLUTION AND SPECTRUM OF THE PERTURBATIONS 

We now proceed to deriving the equation of motion for the gauge- invariant Mukhanov- 
Sasaki variables 

Q 1 = H 1 + ^- (44) 

For the single field, slow-roll case the equations of motions for the variable aQ assume a 
particularly simple form of an oscillator with a time (background only) dependent mass. 
The perturbations in the scalar sector are then easily integrated from an initial state which 
can be approximated by a Minkowski adiabatic vacuum for k ^> aH to a final state at 
horizon crossing k = aH after which the perturbation freezes out and stops evolving. 

For multiple field models two complications can arise. The first is that mixing terms in 
the multiple field potential will lead to a system of coupled oscillators with a non-diagonal 
mass matrix. The second is that a non-trivial field-space metric (kinetic terms) leads to the 
presence of connection terms in the equations. Another important difference with the single 
field case is that, in general, the perturbations will continue evolving once they cross the 
horizon and must be solved for well into the super-horizon regime. 

In what follows we use the spatially flat gauge in which 

^Q = 0, (45) 

a 2 E Q - aB Q = 1*, (46) 

where the subscript Q denotes the use of the spatially flat gauge and \P is the curvature 
perturbation. This gauge makes the equations simple as the Mukhanov-Sasaki variables 
reduce to the simple form: 

Q 1 = H 1 . (47) 



14 



Introducing the choice of gauge into Eq. (|42|) leads the to following expression for Aq: 
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(48) 



Differentiating this we then obtain 
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The value for \I/ comes from Equation (1431) : 
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Thus the curvature perturbation has the expression: 
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(51) 



Using the previously- derived expressions for Aq and $ we can now calculate the gauge 
terms in the right-hand-side of the equations of motion of the perturbation variables, Equa- 
tion flU]). 
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Now we combine all the above results, Equations (j4"9~|) . (I5TI) . and (1521 can now be combined 
to rewrite the equation of motion for each perturbation 5(f) 1 , Eq. ( 133]) . in terms of the new 
variable Q 1 
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H 



r 



ILJ<, 



[KlIlL 



8ttG 



r 



LMY 



G 



IJ<, 



Q J = 0. 



(53) 
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Notice that all coefficients in the equation above depend only on background quantities. For 
a canonical kinetic term, Gjj = Sij, all the connection coefficients vanish and the only mixing 
of the Mukhanov- Sasaki variables Q K comes from the mass matrix of the inflatons. The 
mixing provided by the non-trivial field-space metric enhances the possibility that non-trivial 
features on the power spectrum can be detectable. 

A. Adiabatic and Entropy Perturbations 

Equation ( |53l) allows one to integrate efficiently the system of perturbations from an 
initial state well inside the horizon to a final, super-horizon state at the end of inflation 
when the perturbations determine initial conditions in the constituents of the universe after 
a period of preheating. 

In the single field scenario the perturbations in the inflaton are simply related to per- 
turbations in the curvature. If multiple fields are present however there can exist entropy 
perturbations between the different fields which leave the curvature unperturbed. The two 
types of perturbations induce different initial conditions for the fluid perturbations at the 
beginning of the radiation era and lead to different CMB signatures [261 ] . 

In terms of the generalized trajectory for the variables Q K one can define a single adiabatic 
(or curvature) and n — 1 entropy perturbations at any point along the trajectory. The 
decomposition allows one to identify the observationally distinct perturbations at the end 
of inflation. The variables Q K define a "perturbation vector", Q in field space. We can 
decompose this vector in two components: the one parallel to the velocity of the background 
fields and the component orthogonal to it. We start by defining the unit field velocity vector, 
or the adiabatic direction. 



t = J— . (54) 



The adiabatic perturbation is defined by projecting the perturbation vector Q along the 
adiabatic direction a . 



Gij&Qi ftSft yG k i(p k (j) 1 



Qa = ^—^- = -j== + ^—= V>- (55) 

<G kl ^ a/GWV 
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The corresponding entropy direction is defined as the component of Q orthogonal to a : 



Q - (~a ■ Q)~a 



Q - (~a ■ Q)~a 



The entropy perturbation is the projection of Q along e : 



(56) 



Q t 



Q 



~a ■ Q 



Q - (~a ■ Q)~a 



(57) 



We can now calculate the standard quantities, the co-moving curvature perturbation TZ and 
the re-normalised entropy perturbation, S. 

H 



1Z 



Gijftft 



H 



Qa 
Qs 



(58) 
(59) 



djW 



The corresponding power spectra for the two variables are: 

k 3 



i\(k) 



2vr 2 
2tt 2 



\n{k)\\ 

\S(k)\ 2 . 



(60) 
(61) 



These are the two quantities that we plot for a number of configurations. We look at a 
single brane model, as well as multi-brane models (3 mobile branes) in which the branes 
start either from the same location, or from slightly different positions. We find that the 
spectrum of the adiabatic perturbation is insensitive to the initial conditions and it has the 
same shape for both single-brane and multi-brane configurations. Its shape also agrees with 



the calculations done in Ref. 



11 



211 ] using the Sasaki-Stewart formula [27]. 



VIII. NUMERICAL SOLUTION OF POWER SPECTRA 



To qualitatively study the spectrum of the perturbations produced in these models we 
integrate numerically the evolution of 100 perturbation modes in a model with 3 mobile 
branes. We take an exponential sampling of the momenta such that we have a uniform 
sampling as a function of the number of e-folds throughout the entire inflationary period. 
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The modes are evolved from a time 3 e-folds before horizon crossing. Before that the evolu- 
tion can be solved analytically under the assumption that the mode is evolving deep inside 
the horizon in a deSitter background. The integration is continued until the background 
stops inflating. Fig. [2] shows the evolution of the Mukhanov- Sasaki variable for a single 
wavenumber together with its velocity. The perturbation is initially constant once it crosses 
the horizon but then starts growing due to the roll-over after the inflection point. 
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FIG. 2: Left: Evolution of the Mukhanov-Sasaki variable. The non-trivial shape of the potential 
leads to evolution of the variable after horizon crossing. The late time amplification is due to the 
change in curvature of the potential after the inflection point; the negative mass squared of the 
inflaton leads to a tachyonic amplification of the modes. Right: The corresponding time derivative 
of the same Mukhanov-Sasaki variable. We plot the absolute value of the velocity on a log scale. 
The velocity increases after the inflaton passes the inflection point. To calculate the adiabatic and 
entropy spectra we take a snapshot just before inflation ends. 



A snapshot of the perturbations is taken just before inflation ends and these are then 
decomposed into adiabatic and entropic components following Eqs. (1VII Ah and (IVIIAp to 
obtain the final super-horizon power spectra. 

The single-brane model give the same qualitative results as the multi-brane model in 
which all the brane are initially coincident, and then roll together. We choose the parameters 
such that the adiabatic perturbations are COBE normalized for both the single-brane and 
the multi-brane models. The amplitude of the entropy spectrum is smaller for the single- 
brane case, but the spectra have the same shape. The results are shown in Fig. [3] together 
with the spectral index n s and running dn s /dlnk as a function of e-folds. 



18 



1.5 - 






0.5 " 




0.1 



= 



-0.1 



H 1 H 



H 1 1 1 1 1 h 



H 1 h 



H 1 1 1 1 1 h 





1.5 






- 0.5 



N end -N 



N end -N 



FIG. 3: Upper Left: The adiabatic spectrum has the same shape for single or multi-brane models 
and is insensitive to the initial positions of the branes. Upper Right: The entropy spectrum for 
coincident branes or a single brane. Lower Right: The entropy spectrum for the same model, with 
the branes starting from different locations. The amplitude of the spectrum is still 6 orders of 
magnitude smaller than that of the adiabatic perturbations Lower Left: The spectral index and 
its running for the adiabatic perturbations. The spectrum is blue before the inflaton passes the 
inflection point, so to get a red spectrum as suggested by the WMAP data we need 60 e-folds of 
inflation after passing the inflection point. 



The spectrum of adiabatic perturbations is initially blue before the inflection point and 
then becomes red. The running of the spectral index is small throughout the trajectory. 

The shape of the entropy spectrum is highly sensitive to the trajectory in field space, a 
small change in the initial positions of the mobile branes results in a completely different 
shape of the spectrum as compared with the case of coincident branes. However the ampli- 
tude of the entropy spectrum is 6 orders of magnitude (9 orders of magnitude for coincident 



19 



branes) smaller than the amplitude of the adiabatic perturbations. Such a small contribu- 
tion to an isocurvature component of the primordial perturbations would be unobservable 
be even the most sensitive future CMB experiments. 
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FIG. 4: The marginalized posterior for the spectral index from the "WMAP5yr + CMB" 
combination for the case where running of the spectral index is assumed to be zero. The square 
points are the value of the spectral tilt of the adiabatic perturbations shown in Fig. [3] which crossed 
the horizon at different times during inflation (shown on left axis). Only perturbations crossing 
the horizon at late times are compatible with the red spectra favoured by CMB data. 



In Fig. H] we compare the marginalized posterior distribution for the spectral index ob- 
tained from the "WMAP5yr + CMB" chains 30j (see 28|] for details of the data combination) 
with the spectral index of the perturbations which crossed the horizon at different e-foldings 
from the end of inflation (iV on d). The chains were obtained for a pure power law model 
of adiabatic primordial perturbations. Only the modes which exited the horizon after the 
inflection point are compatible with the red spectrum preferred by the data. 
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The comparison for when running of the spectral, n mn = dn s /d\nk index around a pivot 
point is allowed is shown in Fig. [5] for the same data combination. The pivot is at k = 0.002 
Mpc -1 scales probed by current CMB experiments and the marginalized confidence contours 
are shown for the n s , n run plane. The smallness of the running observed in the models is in 
slight disagreement with the data which prefers a negative running of the index at around 
the 2-cr level. The shift in the preferred n s value to n s p» 1.05 is due to the scale dependence 
of the spectral tilt in the presence of non-zero curvature. 
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FIG. 5: The adiabatic perturbations results of Fig. [3] plotted over the CMB constraints when 
both n s and dn s /dlnk are fitted. The 68% and 95% confidence regions are obtained from the 
"WMAP5yr + CMB" 28J. The blue squares on the trajectory indicated, from left to right, the 
point at 20, 30, 40, 50, 60, and 70 e-foldings from the end of inflation. The smallness of the spectral 
running in the multi-brane models leads to some tension with the observations which prefer more 
negative running. 
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IX. CONCLUSIONS 

We analysed the possible signatures that brane-inflation models could give, that may be 
detectable in observations of the CMB spectrum. To do so we derived the evolution equation 
for the perturbations of the inflatons and metric and numerically evolved a large number of 
modes starting inside the horizon, all the way to the end of inflation. We then separated the 
fluctuations into adiabatic and entropy ones and derived their respective spectra. We found 
that the adiabatic spectrum is insensitive to the inflaton trajectories and has the has the 



same shape as was determined in earlier work III, |21[ using the Sasaki- Stewart formalism 



27| . The entropy one is highly sensitive to the inflaton trajectories but its amplitude is too 
small to be observable. We also offered an analytical argument that in a large class of models 
where most of the e-folds come from a small region around an inflection point of the potential, 
the inflaton trajectory is almost a straight line and therefore the entropy perturbations are 
expected to be small. Thus multi-brane models solve a number of fine-tuning problems but 
due to their e-foldings coming from a reduced section of the field trajectory they display less 
phenomenology than single field models. This argument would lead us to expect that other 
potentially observable effects such as non-Gaussianities would also have a low amplitude 
in these models. Typically large values of the non-Gaussian / NL parameter are produced 
during periods of high acceleration in the field trajectory which are not present in these 
models 10j. 

Interestingly, the old models of brane inflation using branes oriented at angles and toroidal 
compactifications |2j, y|, U, |29] have highly non-trivial inflaton trajectories as a result of the 



brane interaction inside a compact manifold. This leaves open the possibility that more 



complicated models, like multi-throat models 23] will feature non-trivial inflaton trajectories 



and therefore observable features in the power spectrum and non-Gaussianity. 

As we have shown the models have no particular problem in reproducing the observed 
tilt in the power spectrum but as expected in models that probe short trajectories, the 
running is small. If the observed running is confirmed to be at the current best-fit value 
of dn s /d\nk ~ —0.05 with increasing significance in future measurements it may turn out 
to be the strongest constraint on these models. As with most multi-field inflation models a 
detection of a gravitational wave background would also be in conflict with the models since 
the tensor perturbation are not boosted by the number of fields contributing to inflating the 
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universe. 
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APPENDIX A: DETAILS OF THE CALCULATION 
1. The F-term potential 

The F-term is calculated with the usual formula: 



V F = e K (K ab D a WD b W - 3 \W\ 2 ) , 



(Al) 



where 

/C = -31og(r) 

The covariant derivatives are expressed as: 



(A2) 



D T W = d T W + Wd T lC = d T W - -W 

r 

df 3 



D T W = chW + WchJC = djW 



-W 



Qt. 3 

D&W = chrW + Wch-K, = -£--W 



>T r 



(A3) 
(A4) 
(A5) 
(A6) 



The inverse Kahler metric needed in the expression above, can be calculated analytically in 
the general case: 



/ j. I dfl dfl I d 2 fl I I df n df n I d 2 f n dfl / d 2 fl 9/2 / d 2 f 2 



/c 



ah 



dfl I d 2 h 
d<fii I dipidifi'^ 

dtp2 I d<fi2d<P2 



9fn j d 2 f n 



1/ 92 h 

' dipid(fil 

l/J^- 



dip2dip~2 



df n j d 2 f n \ 






1/ 92 fn , 

' d<p n d<fi^ / 
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Using the above expression for the inverse Kahler metric one obtains the following terms: 

dfn dfn \ 

\ d<t>id<t>l 

/ au 

3, 






h' Ti ■■Tmvd \v = - I chW -~w)w (-— m +■■■ + -— aK 

UTWU^W \OtVV wjw I ^ d2fi i- i" ^ _^_ 

V / \ T' 1 3J,. flit ™" «,a ^^,a 



v / \ vi a0 l9 0* vn ^ n9( 

K^DTWD, W = -WW (l^M _J_ + . . . + i^Mi 4_ 



(A8) 

(A9) 

(A10) 



Adding all the above terms and —3 |W| we obtain the F-term: 

dh-dfl. dfn djn. \ 

<h W ( ) T Y\ • J .1 ■ - -^M- + • ■ ■ + °+^- - r (W^PF + Wd T W) 



V, 



oio<p 1 



(All) 



We assumed that the superpotential itself is independent of the fields 0i . . . <p n . 
For a racetrack superpotential (jHJ) the corresponding F-term potential becomes: 



V, 



aA 2 e~ acj a 



3r 2 <92 ^_ d 2 3_ 
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abAB (^lj± d _i + ^X^^£ _ r d 2 f_ d 2 g \ 
\dipdip dcp dcp dcpdcp dip dip T 'dcpdcp 1 dipdip ) 



dipdip 
cos ((a — b) u) 



(A12) 



The expression above becomes much simpler if the functions / and g are taken to be the 
simplest ones, / = cpcp and g = ipip. 

2. Perturbation of the kinetic term 



This is the first term in Eq. ( |25l) : 



-A W=ggTQu<iF\ 



(A13) 
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and we expand this term as follows: 
1 



--d„ 



det (g + 8g) (<T + Sg^) d v (0 K + 5cj> K ) 



K i xaK} 

, 1 \ / I I i 'I . 1 / / — t ) t I I [ if Oil' I ( f . . I ( : ) ' 

y/- det (g + 6g) 
=d» y^~gg^dA K \ +-^=d ll y=g-g^dM K ] + ~^d, [y/^gSgT 'd v <j> K ] 



background 



-.d„ 



^Tg-Tr [g-%] g^d, 



>" v %(j) K 



1 1 



-9 2 



-Tr [g- l 5g] d, [^ g ^d v cj) K ] (A14) 



In the above equation we used the expansion of a determinant: 

det [g + Sg] = det [g (l + g~ x Sg)] = det g det [(l + g'Hg)] = detge Trlog ( 1+g ~ 159 ) ~ 
det ge Tr \. 9 ~ Ha \ ~ det g(l + Tr [g- l 5g] ) (A15) 

and therefore: 

det [g + 5g] = ^-det^ f 1 + ^Tr [^fy] J . 



(A16) 



We now proceed to simplify each one of the 4 terms in Eq. (jA14jl above using the back- 
ground metric Eq.( l22|) and the fact that the background fields and metric depend only 
on time. We write the Laplacian of a generic field $ in terms of the momentum as 
d m d m & = — fc 2 $ and proceed with the details of the perturbations of each of the 4 terms in 
Eq.dXTU): 

a. 1st term 



This term is straightforward to compute, as it involves only the perturbation of the scalar 
field: 

\/—g a 6 



■^d 
a 6 



3 A^k 



a A 5(t) 



+ -,d m d m 86 K 



a a 2 



(A17) 



b. 2nd term 

In order to proceed we will need the form of the perturbations of the inverse metric. We 
can easily find the inverse metric imposing the condition that 



(9, 



\iv i u yfj,uj 



9 vX + 6g vX ) = 6* 



(A18) 
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to first order in the perturbations. We find the inverse metric to have the expression: 

-(l-2A(i,x)) J(t)B,i(«,x) 



g^ + 5g 

' \ (t) B, (t, x) 4, (i) [(1 + 2^ (t, x)) ^ - 2£„ (t, x)] . 

With this expression the 2nd kinetic term perturbation becomes: 
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a a 



c. 3rd term 



This terms has the expression: 
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^-g-Tr [g-Hg] g^d A 



.K 
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-(j) K -Tr [g l 5g\ - <p K -Tr [g l 5g] - —<p K Tr [g 1 8g] 

where the trace of the metric perturbation has the form:: 
1, 



a 3 -TV [g x 5g] 



Sii 



-Tr [g- l 5g] = A - 3^ + B u + £ 22 + E, 33 = A - 3^ - k 2 E. 



(A21) 



(A22) 



d. 4th term 



Finally, this is the simplest of the 4 terms and it has the expression: 



1 }-Tr [g-'Sg] d, [v^jT^*] 



U'KTr [g-Hg] + ^Tr [g'Hg] 



2a 3 



Tr [g-Hg] d 



a 3 (-1)0 



A" 



(A23) 



2 r L " " J 2a 

We see that these two terms we obtained will cancel 2 of the terms from the previous 
expression. We can now move on to the potential terms. 



3. Details of the Calculations: M-S variables 



We present here the details of the calculation for the term: 



9A r KidV 
~ 2A Q G 7u7 + 



A lk \ 



hK 



(A24) 



26 



appearing on the R.H.S. of the perturbation equation. We want it to take the appropriate 
form in terms of the Sasaki- Mukhanov variables. We will therefore use Eq. (j4"9l) and fl5Tl) in 
the expression above as well as the equation of motion for the background fields, Eq. (1231) . 
The expression above becomes: 
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(A25) 



We use the equation of motion to express J^ in terms of the fields and their derivatives. 
We also notice that the terms Gu<p I 6<p J cancel and the terms containing the derivatives of 
the Kahler metric can be combined to form the Christoffel symbol with all indices lowered. 
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(A26) 



Finally, the term ^-G LM 4 >L <i) M G I j4> I 5(j) J can be related to H/H. We start out by calcu- 
lating H/H: 
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The above result tells us that we can substitute: 
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The expression we want to calculate becomes: 
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In the end we obtain the equation: 
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